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ABSTRACT 

We present the highest redshift detections of resolved Lya emission, using Hubble Space Telescope/ACS 
F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science pro- 
gram in the GOODS CDF-S. We detect Lya emission from three spectroscopically confirmed z = 4.4 Lya 
emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lya emission. Comparing 
the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the 
escape of Lya photons at high redshift. While our data do not support a positional offset between the Lya and 
rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are signif- 
icantly larger in Lya than in the rest-frame UV continuum. This result is confirmed when comparing object 
sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is signifi- 
cantly less than observed in similar filters from the ground. These results together imply that the Lya emission 
is not strictly confined to its indigenous star-forming regions. Rather, the Lya emission is more extended, with 
the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lya 
photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing 
a significant role in galaxies at these redshifts. 

Subject headings: galaxies: high-redshift - galaxies: evolution - galaxies: ISM 
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L INTRODUCTION 

High-redshift Lya emitting galaxies (LAEs) are some of the 
most intriguing objects in the distant universe. Their strong 
Lya emission was thought to be indicative of the first galaxies 
(iPartridge & PeeblesI 19671) . implying that they could possibly 
contain the first stars and likely be composed of pristine gas. 
However, recent studies of their physical properties imply 
that some of these galaxies may be more evolved, with many 
LAEs exhibiting rest-frame ultraviolet colors indicative of 
modest-to-moderate dust extinction (e.g.. 'Pirzkal et al .ll2007l: 
Lai etini2007l: iF inkelstein et al.ll2008. 2009; PentericcTSian 
2009HOno etSLll2010 i. As Lya photons are resonantly scat- 
tered by neutral hydrogen, galaxies with dust would be un- 
likely to exhibit Lya in emission. Thus, just how Lya escapes 
from a galaxy with a dusty interstellar medium is an outstand- 
ing question in the study of distant galaxies. 

As LAEs have been selected on the basis of their Lya emis- 
sion, some mechanism must allow the escape of these pho- 
tons. One possibility is that the Lya photons we see have been 
shifted out of resonance by scattering off of the receding edge 
of an outflow in the ISM. Evidence for outflows has been ob- 
served many times in the t ypically rnore evolved Lyman break 
galaxies (LBGs; e.g., Sha plev et alj|2003l: iBielbv et al.ll2010t) 
as a velocity difference between Lya emission and ISM ab- 
sorption features. O nly recently hav e outflows been shown 
to exist in LAEs, as iMcLinden et al] (|2pi0) discovered that 
in two LAEs at z ~ 3.1 Lya emission had a slightly higher 
redshift than the rest-frame optical [O III] emission, which is 

^' Cerro Tololo Inter- American Observatory, La Serena, Chile 
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thought to come from H II regions at the systemic redshift. 
In either case, much of the Lya emission is shifted redward 
of the resonance line at 1216 A, and thus will have an easier 
chance of escaping, even in a uniform ISM. 

Alternatively, if the line emission is primarily at reso- 
nance, much of the Lya emission can still escape if the 
ISM is primarily clumpy, as the Lya photons will scat- 
ter off of the clumps, and b e screened from seeing much 
of the dust (iNeufeldl [1991 ilans_en & Oh 2006). This 
type of ISM geome try can explain the dustiness of LAEs 
at z 4.5 (Finkelstein et alj 120081 |2009) . and can also 
explain the large number of high L ya equivalent widths 
(EWs) which have b een observed ( e.g..lKudritzki et ^1120001 : 
iMalhotra & Rhoads|[2 002: Finkelst einet al.ll2007h . 

In either of these radiative-transfer scenarios, any detected 
Lya emission will be spatially de-correlated from its origina- 
tion point within its host galaxy. By comparing Lya emission 
from a narrowband filter to the rest-frame UV emission from a 
neighboring broadband filter, one can diagnose whether this is 
the case; if Lya has undergone any extreme radiative-transfer 
effects, this should reveal itself in a larger size in the Lya 
emission, as well as possibly a diffuse Lya halo. 

However, all currently known LAEs have been discovered 
via ground-based narrowband imaging, which even in the best 
seeing conditions cannot resolve the extremely small physical 
sizes of LAEs of 1 - 2 kpc (Bond et al. 2009, Malhotra, S. 
et al. 2010, in prep). Here we report on the results of a new 
Hubble Space Telescope imaging survey designed to search 
for resolved Lya emission from LAEs at z 4.4 using the 
F658N narrowband filter on the Advanced Camera for Sur- 
veys. Throughout we use the AB magnitude system, where 
m^B = -2.5 log (f^) - 48.6 mag. Where applicable, we as- 
sume a concordance cosmology, with Ho = 70 km s"' Mpc"' , 
flm = 0.3 and = 0.7. At z = 4.4, this corresponds to an 
angular scale of 6.671 kpc arcsec"^. 

2. DATA 

2.1. Observations 

Thanks to the successful repair of the Advanced Cam- 
era for Surveys (ACS) during Servicing Mission 4 (SM4) 
to the Hubble Space Telescope (HST), we were able to ob- 
tain ACS parallel imaging during the Early Release Science 
(ERS; Win dhorst et al.l i2010) Wide Field Camera 3 (WFC3) 
observations of the Great Observatories Origins Deep Survey 
(GOODS) Chandra Deep Fields - South (CDF - S). We ob- 
tained 1 1 orbits per pointing over 8 independent pointings. 
Due to the location of the WFC3 fields, all of the ACS point- 
ings overlapped the GOODS CDF-S field, which has existing 
deep public data in the F435W, F606W, F775W and F850LP 
ACS filters (as well as a wealth of other multiwavelength 
data). A detailed summary of the primary W FC3 ERS images, 
as wel l as their layout and analysis is given by lWindhorst et all 
(l2OT0h . 

We split each parallel pointing into 9 orbits with the F658N 
narrowband filter and 2 orbits with the F814W broadband fil- 
ter. With a central wavelength of 6584 A and a full- width 
at half-maximum (FWHM) of 73 A, the F658N observations 
will detect Lya (which has X„s, = 1215.67 A) from redshifts 
4.386 < z < 4.445. A model spectrum of a LAE at z = 
4.42 is shown in Figure [T] At z = 4.4, the existing GOODS 
ACS data cover rest-frame wavelengths of ^ 800 A (F435 W), 
1100 A (F606W), 1400 A (F775W) and 1600 A (F850LP). 



0.5 
0.4 

I- 0.3 

e 0.2 

0.1 
0.0 











-i 






1 


\ 


I.... 1...: 



4000 6000 8000 10000 
Lambda (A) 



Fig. 1 .— A model spectrum of a z = 4.42 LAE, showing the HST ACS 
F435W, F606W, F775W and F814W bandpasses in blue, green, red and 
brown, as well as the F658N bandpass (dotted Hne). From z = 4.38 - 4.45, 
LyQ passes through the F658N bandpass, allowing imaging of Lya light at 
high redshift. 

The F814W data will cover the continuum at ^ 1500 A, pro- 
viding an independent observation in addition to the existing 
GOODS dataset. 

2.2. Data Reduction 

The raw ACS data were downloaded from the Space Tele- 
scope Science Institute (STScI) archive. The ACS data were 
taken in 31 separate visits, with typical dithers within each 
visit of < 20". Images from visits with central pointings sep- 
arated by less than 25" were reduced together, yielding 17 
separate reductions. The raw images were processed using the 
calacs task, which is in the stsdas package in IRAF'. This 
task provides routine calibration, including bias, dark and flat- 
field corrections, using the most recent ACS reference files 
taken after SM4 retrieved from the HST archive. ACS data 
obtained after SM4 suffer a low-level striping pattern. We im- 
plemented a custom-built script (provided by NAG) to remove 
this pattern prior to the flat-field correction. 

The calibrated and pattern-corrected images were cleaned 
of cosmic-rays, distortion-corrected , registered and co mbined 
using the task multidrizzle dKoekemoer et al.l 120021) . 
Upon completion of the initial run of multidrizzle in 
each visit, it was apparent that the registration was not ideal, 
as stars in the combined images appeared elongated. We thus 
ran custom-built scripts (provided by AMK) on a visit-by- 
visit basis to correct the World Coordinate System (WCS) in 
the headers of the individual frames, solving for the relativ e 
astrometric shifts between frames dWindhorst et alj l2010h . 
Multidrizzle was then run a second time to create a fi- 
nal, combined image for each visit. As we planned to use 
the existing GOODS ACS data in our analysis, we used the 
GOODS ACS image sections as reference images when run- 
ning multidrizzle. In order to correct for small astro- 
metric differences between the GOODS data and these new 
ACS data, we first ran multidrizzle with no reference 
image, and performed photometry using t he source extrac- 
tor software package (hereafter SExtractor: lBertin & ArnoutsI 
61996) ) to identify objects in the image. We did the same 
to the relevant GOODS section, and then ran the IRAF tasks 
xyxymatch and geomap to match common objects be- 

' IRAF is distributed by the National Optical Astronomy Observatory 
(NOAO), which is operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation. 
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Fig. 2. — F658N (left) and F814W (right) image of GOODS-S section 23. This was the only GOODS-S section which was completely covered by our 
observations. The narrow bandpass of F658N is apparent when comparing the depths of these images, as many fewer objects are apparent in F658N even though 
the exposure times were longer. 



tween the two frames, and construct a shift file. This shift file 
was used for the final iteration of multidrizzle to create 
F658N and F814W images matched to each GOODS section 
covered by our observations. Our final dataset was composed 
of one image in each of the two filters for the 1 1 GOODS-S 
sections that we covered: 12, 13, 14, 22, 23, 24, 32, 33, 34, 
42 and 43. Typical exposure times in the reduced datasets are 
11000 s in F658N and 2200-2500 s in F8 14W. Images of sec- 
tion 23 in the F658N and F8 14W bands are shown in Figure|2l 

2.3. Catalog Construction 

We created narrowband-selected catalogs for each observed 
GOODS-S section using SExtractor in two image mode, with 
the F658N image for each section as the detection image, 
and our F658N and F814W images, as well as the GOODS 
F435W, F606W, F775W and F850LP images as the mea- 
surement images. We used identical SExtractror parame- 
ters as used in GOODS. The final catalog encompassing all 
covered sections includes 3081 narrowband-selected objects, 
with fluxes measured in 0.7" diameter apertures, as well as 
estimates of the total flux using SExti-actor's MAG_AUTO 
measurement. 

Initial flux errors were taken to be the calculated SExtractor 
errors. We checked these errors by measuring our own errors 
in each image. This was done by measuring the flux in lO** 
randomly placed 0.7"-diameter apertures in each of the six 
images, and then examining the spread of these fluxes (this 
was done in GOODS-S section 23, as this was the only sec- 
tion which had complete coverage by our F658N and F814W 
data due to the unfavorable positioning of the parallel expo- 
sures with respect to the GOODS sections). The characteristic 
1 a error for each image was taken as the cr of a Gaussian fit 
to a histogram of the flux distribution. Comparing this error 
to the median SExtractor error in each image, we find that 
SExtractor underestimated the errors by up to ^ 20% (with 
the exception of the F814W data, where SExtractor overesti- 
mated the errors by 18%). While we trust our independently 
computed errors as being indicative of the global uncertainty 
in the image, the errors computed by SExtractor include infor- 
mation on the local background. We thus scaled the median 
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Fig. 3. — Number counts of objects in the F658N data. The colored his- 
tograms denote the numbers in each section, with the values being given by 
the left-hand vertical axis. This thick black line denotes the total combined 
number counts, with the values being given by the right-hand vertical axis. 
The varying number of objects per section is related to the amount of section 
area which received F658N coverage. The dotted line denotes where the total 
number counts fall to 50% of their peak value, which is at m^ sjgjv = 25.3. 

SExtractor error to match the global uncertainty in each band. 
The derived 5 a limits for each band in a 0.7"-diameter aper- 
ture are: 25.0 (F658N), 27.1 (F814W), 27.4 (F435W), 27.6 
(F606W), 27.0 (F775W) and 26.8 (F850LP). Number counts 
of objects in the F658N images are shown in Figure[3] show- 
ing a peak at mfesgw ^ 24.8 mag. 

3. SAMPLE SELECTION 

3.1. Spectroscopically Confirmed LAEs 

In lFinkelstein et al.1 (120081) and iFinkelstein et al.1 (l2009t) . a 
sample of 14 LAEs were discovered in the GOODS-S field 
using ground-based narrowband selection. These studies used 
three overlapping narrowband filters, centered at 6560 (here- 
after NB656), 6650 (NB665) and 6730 (NB673) A to discover 
LAEs at z « 4.4 - 4.5. Samples of 4, 2 and 8 candidate LAEs 
were discovered in the three images, with 5 cr depths of 24.9, 
25.0 and 25.2 mag, respectively. The ACS F658N filter can 
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Fig. 4. — Cutouts of the three LAEs in our sample, 2" on each side. The black circles are centered on the F658N centroid, with a 0.7" diameter The F658N 
data is from our program, while the remaining images are from the GOODS dataset. Object CHa-2 is near an image edge, thus the lower-right comer of this stamp 
is not real data. Although trailing due to charge transfer inefficiency is apparent in object CHa-2 (as blurring from lower-left to upper-right), the narrowband flux 
is still detected at 3.9 cr in a 0.7" diameter aperture. The second column shows the CTE-corrected stamps, as discussed in §5.2.4. The CTE-correction needs to 
be verified with newer calibration data before it can be folded into the analysis, but improvement in the background can be seen, especially near LAE CHa-2. 
The F435W data probes rest-frame A ~ 800 A at z =4.4, thus the non-detections in this image are expected. The cyan circle in the F606W images denotes the 
hya emission centroid after connection for a systematic positional offset between the WCS of the new F658N data and the existing GOODS data. The red circle 
denotes the 1 a uncertainty on the positional offset. These uncertainties are large in CHa-2 and CH8-2. However, the offset can be computed to a high precision 
in CH8- 1 , and we find that the apparent offset between the hya and UV emission can be explained by WCS differences in the two datasets. 



TABLE 1 
Properties of the LAEs 



Object 


Redshift 


RA 
(J2000) 


Dec 
(J2000) 


(mag) 


(mag) 


m;r775iv 
(mag) 


^PER 
'*,f658iV 
(kpc) 


,APER 
'/i.f775lV 
(kpc) 


Rest EW 
(A) 


CHa-2 
CH8-1 
CH8-2 


4.414 
4.434 
4.433 


03:32:39.77 
03:32:49.01 
03:32:54.04 


-27:51:14.97 
-27:49:02.08 
-27:50:00.83 


24.90 ± 0.28 
25.20 ± 0.22 
24.64 ± 0.20 


26.57 ± 0.09 
27.06 ±0.14 
25.93 ± 0.05 


25.56 ± 0.06 
26.54 ±0.14 
24.98 ± 0.04 


0.73 ±0.13 
1.40 ±0.16 
1.33 ±0.18 


0.67 ± 0.05 
1.07 ±0.11 
0.67 ± 0.03 


167 
176 
53 



Note. — The position is the centroid of the F658N counterpart. All magnitudes were measured with 0.70" diameter apertures. Half-light 
radii (r;,) were converted from arcsec to kpc assuming all objects are a t z = 4.4^ which gives an angular scale of 6.671 kpc arcsec"' for our assumed 
cosmology. The rest-frame equivalent widths are from , Finkelstein et aUf20()9l) . and are based on the ground-based narrowband imaging. 



measure Lya emission from galaxies at z = 4.38 - 4.45; thus 
it would also observe Lya from objects discovered in the red 
half of the NB656 filter, or in the blue half of the NB665 fil- 
ter Of the six candidate LAEs discovered in these two fil- 
ters, three fall in the area covered by our F658N observations. 
These three objects are CHa-2 , CH8 - 1 and CH8-2, using the 
nomenclature from Finkelstein et alj (l2009h (where CHa de- 
notes CDFS Ha, i.e., NB656, and CH8 denotes CDFS HaH-80 
A, i.e. NB665). Images of these three objects are shown in 
Figure |4] 

Although these objects were previously selected via nar- 
rowband observations from the ground, they had yet to be 
spectroscopically confirmed. We recently obtained optical 



spectroscopy of these three objects with the Inamori Magel- 
lan Areal Camera and Spectrograph (IMACS) at the Magel- 
lan Baade Telescope on 11-12 November 2009 (NOAO PID 
2009B-0371, PI Finkelstein). The full details of this spec- 
troscopic dataset will be presented in a future paper (Zheng 
et al. in prep), but in brief, each object was observed as part 
of a 4-hour slit-mask integration, with the f/2 camera and the 
300 lines/mm grating blazed at 17.5° (giving R ~ 1000). The 
reduced, one-dimensional spectra of these three objects are 
shown in Figure|5] Each object exhibits a single emission line 
with no significant continuum light, indicative of Lya emis- 
sion at high redshift. Fitting a Gaussian curve to these emis- 
sion lines, we find redshifts of CHa-2, CH8-1 and CH8-2 of 
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4.414, 4.434 and 4.433, respectively, placing the Lya emis- 
sion line of each object in the bandpass of the ACS F658N 
filter. 

Examining these objects in Figure ID they all appear ro- 
bustly detected in the F658N image. However, especially in 
the case of CHa-2, the noise due to poor charge transfer effi- 
ciency (CTE) in the nearly decade-old CCDs onboard ACS is 
apparent. Nonetheless, when we consult our narrowband se- 
lected catalog, we find that CHa-2, CH8-1 and CH8-2 are all 
formally detected, with detection significances of 3.9, 5.0 and 
5.4 (T, respectively. Combined with the fact that LAEs were 
previously known to reside at these locations, we are confi- 
dent that we are in fact detecting Lya emission with ACS. 

3.2. Photometric Redshift Selection 

In addition to objects previously selected on the basis 
of their Lya emission at z « 4.4, we have also exam- 
ined the F658N images for objects which are likely to re- 
side at z ^ 4.4 based on their spectral energy distribu- 
tion. We selected objects at this redshift from two cata- 
logs, both from S. Cohen et al. (2010, in prep). The first 
consists of ^ 15000 objects with spectro-photometric red- 
shifts computed using both ACS broadband and grism slitless 
spectroscopic data from the Probing Evolution and Reion- 
ization Spectroscopically (PEARS) program (PI S. Malho- 
tra). The second catalog consists of ^ 8000 photometric 
redshifts measur ed over the e ntire GOODS-S region, using 
VLT/VIMOS U (iNonino et al.1.2009) . GOODS/ACS v2.0 B, 
y, i', z', and GOODS VLT/ISAAC v2.0 J, H and K-band data 
jRetzlaffet al.ll20Tol) . 

We examined these objects for galaxies with best-fit 
(spectro-)photometric redshifts of 4.38 < z < 4.45, placing 
any Lya emission in the F658N bandpass. We also included 
objects that had this redshift slice contained within the 68% 
confidence range on their spectro-photometric redshift. We 
found 106 objects meeting these criteria. We then matched 
these objects to our F658N catalog, using a matching radius 
of 0.5", and we found six objects that have F658N counter- 
parts. The low number of matched objects is expected, as 
only galaxies exhibiting Lya emission at the specific redshift 
placing it in the F658N bandpass would be detected in the 
narrowband data. These objects were visually inspected in 
the F658N data. Of these six objects, only two have mod- 
erate narrowband excesses {vapmew - mfessw = 0.6 and 1.8 
mag). However, both objects have significant detections in 
the F435W-band. This band is entirely blueward of both the 
Lya and Lyman continuum break at z « 4.4 (see Figure [T]); 
thus there should not be a F435W detection if these objects 
were truly at z « 4.4. We conclude that these two objects are 
low-redshift interlopers, and we exclude them from further 
study. 

4. RESULTS 

With our sample of three F658N-detected high-redshift 
LAEs, we investigate their light profiles, as well as the lo- 
cation of their Lya emission. 

4.1. Positional Differences Between Lya and Rest-Frame UV 

Emission 

In a number of objects the Lya emission appears offset 
from the centroid of the rest-frame UV emission. If this ef- 
fect is real, it is quite interesting, as it could indicate that Lya 
is escaping only after scattering off of gas or dust outside the 
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Fig. 5. — Magellan/IMACS optical spectra of the three LAEs in our data, 
centered around the observed Lya emission line. The vertical scale is in 
arbitrary flux units. From the position of the emission line, and the lack of 
any other emission line in the full spectra, we confirm that all thi'ee objects 
are LAEs at z 4.4. This confirms that the flux we detect in the F658N 
image is Lya emission from these objects. 

prim ary stellar populat ion, per haps due to outflows in the ISM 
(e.g.. lWindhorst et al ." 1998: Waddingtonetan[T99l . How- 
ever, we first need to investigate if the offset is real, or if it is 
an artifact of mismatches between the WCS of the new F658N 
data and the existing GOODS-S data. 

We investigated these offsets by examining the relative 
pixel positions of all objects in the images around the LAEs. 
To find these objects, we first ran SExtractor on both the 
F658N and F606W images, using each image as its own de- 
tection image, such that we obtained object coordinates native 
to each image. On a LAE by LAE basis, we first searched the 
F658N catalog for all objects in a given section, excluding 
objects near the edge of our images, as well as objects be- 
low the point where the number counts fall to 50% of their 
peak value, which is at 25.3 mag. We then computed the dis- 
tance in pixels from the LAE to each of these objects. We 
selected objects within a threshold radius, which ranged from 
500 - 3000 pixels in 100 pixel increments, and matched them 
to objects in the F606W catalog, keeping objects that were 
matched within 20 pixels (which is larger than the largest ap- 
parent shift; see FigureH]). By including only objects near the 
LAE, we ensure that we are locally measuring any offset be- 
tween the F658N and F606W image frames. At a 1000 pixel 
radius, on average a dozen matches were found, increasing 
to ^ 40 matches by 2000 pixels. The pixel offsets were then 
computed as the mean difference between the narrowband po- 
sition and the broadband position for each of the matched ob- 
jects. An estimate of the uncertainty on these shifts was taken 
to be the standard deviation of the positional differences for 
the matched objects. 

Figure|6]shows an example of this process, showing the re- 
sults for LAE CH8-1. We plot lines showing both the pixel 
offsets, as well as the offset uncertainties as a function of 
search radius. We chose offset values for each object to be the 
pixel offset value at the radius where the offset uncertainty 
was a minimum. For this object, the x-offset is negligible 
(0.24 ± 0.98 pixels), while the y-offset is significant, at 3.25 
± 0.73 pixels. Investigating FigureHl this offset would move 
the broadband counterpart down vertically, bringing it more 
in line with the narrowband position. This is shown by the 
cyan circle in Figure |4] A similar process was done for the 
remaining objects, and these offsets are tabulated in Table 2. 
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Fig. 6. — The derived pixel offsets between the F658N data and the 
GOODS F606W data for objects near LAE CH8-1, versus the search radius 
used to find objects to derive the offset. Red and blue denote the x and y off- 
sets, respectively, while the solid and dashed lines denote the offset and offset 
uncertainty. Offsets were derived using objects within 500 - 3000 pixels, in 
100 pixel increments. The value of the offset was defined to be the offset de- 
rived from the search radius which produced the smallest offset uncertainty 
(designated by the dotted lines). In this object, the offset uncertainties for 
both x and y reached a minimum at 900 pixels, thus the pixel offsets were 
taken to be the values at that search radius. At larger search radii, the offset 
uncertainties increase dramatically. 

In Figure |4] we show the corrected Lya emission position 
by a small cyan circle, and the typical offset error as a red 
circle (where the radius of the circle is the mean of the x- 
and y-offset errors for a given object). We find that in CHa-2 
and CH8-2, the uncertainties on the derived pixel offsets are 
large, and thus any apparent offset would be at low signifi- 
cance. However, in CH8-1, the offset uncertainties are small, 
and we can see that the computed offset is consistent with the 
centroid of the UV emission. Thus, while the Lya and UV 
emission in CH8-1 appear to be offset, this can be explained 
by relative offsets between the two datasets. 

We conclude that while it is possible that ISM scattering 
effects can result in an offset between the apparent positions 
of Lya and rest-frame UV emission, we cannot conclusively 
support this with our data. The most convincing offset is in 
CH8-1, as the offset is large. However, this is also the only 
object where the offset uncertainties are small enough that we 
can reasonably correct the Lya emission position, and we find 
that its corrected position is then coincident with the U V emis- 
sion. In addition, when inspecting the F814W data taken at 
the same time as the F658N data, the apparent offset between 
F658N and F814W is much less. We move forward assuming 
that the Lya emission is coincident with the rest-frame UV 
emission in all objects. 

4.2. Physical Size of Lya Emission 
4.2. L Individual Objects 

In order to measure the physical sizes of the LAEs in our 
sample in both their Lya and rest-frame UV co ntinuum light, 
we employed the method of iBond et al.l (120091) . We first cut 
out 101 X 101 pixel (3"x3") postage stamp FITS images cen- 
tered on each LAE in both the F658N and F775W data (we 
used the F775W rather than the F606W data for the rest-frame 
UV as it is completely redward of Lya at this redshift). We 
then ran SExtractor on each stamp, using the stamp as both 
the detection and measurement image to determine the flux- 
weighted center of the object, as well as to determine whether 



TABLE 2 
Corrections to Narrowband 
Emission Position 



Object 


Ax 


Ay 




(pixels) 


(pixels) 


CHa-2 


-1.17 ±5.10 


-0.08 ± 4.66 


CH8-1 


0.24 ± 0.98 


3.25 ± 0.73 


CH8-2 


2.74 ±4.16 


2.09 ± 2.28 



Note. — The derived pixel correc- 
tions to the F658N emission position due to 
WCS differences between the F658N and 
the GOODS data. The corrected narrow- 
band emission position is given by the cyan 
circles in Figure|4] while the positional un- 
certainties are shown by the red circles.. 

an object is made up of sub-clumps that might have been split 
up by SExtractor. Experimenting with various values of the 
DEBLEND_NTHRESH parameter, we found that none of our 
LAEs can be split into multiple objects; thus we conclude 
that all objects in our sample are composed of single domi- 
nant components. Inspecting the SExtractor results, we find 
that (as expected) each object is detected in both F658N and 
F775W. 

Using the SExtractor-derived center in the F658N and 
F775W images, respectively, we measured the flux in a series 
of 32 apertures, with radii ranging from 0.015 - 1.2" using 
SExtractor. In order to ensure that the CTE-affected back- 
ground was subtracted as well as possible, we manually sub- 
tracted the background prior to running SExtractor, using the 
iterative mean computed with the IDL task djs_iterstat.pro. 
We then forced SExtractor to assume a background value of 
zero. Previously measured half-light diameters of LAEs are 
- 0.2 - 0.4" Bond et al. (2009, 2010), thus we assume that 
the flux at a radius of 0.6" approximates the total flux. We 
then compute the radius at which the flux is half of the flux 
at r = 0.6", and use that as an estimate of the half-light radius 
(r/,). 

Values of r/, were computed for each object in F658N and 
in F775W, and are tabulated in Table 1 . The curves-of-growth 
(CoGs) of each object in both bands are shown in Figure |2l 
We also show the uncertainty in the CoGs as the shaded re- 
gion, using a similar exercise as explained in §2.3 to com- 
pute the flux uncertainty in all 32 apertures (using the flux un- 
certainty from the appropriate image section). We then used 
these errors to compute uncertainties on our derived half-light 
radii by running a series of lO"* Monte Carlo simulations. In 
each Monte Carlo simulation, we vary the flux at each point in 
the CoG by a random number (drawn from a Gaussian distri- 
bution centered at zero with a = I) multiplied by the flux un- 
certainty, and rederive the half-light radius. The uncertainty 
on the radius is then the standard deviation of the radii from 
the simulations. The radii uncertainties are shown as cross- 
hatched regions in Figure|2] Typical uncertainties on r/, are 
0.02" in F658N and - 0.01" in F775W. 

To determine whether a particular object is resolved, we 
performed the above analysis on a point-spread-function 
(PSF) made from stars in both the F658N and F775W data. 
These PSFs were made by adding together images of five 
stars identified in section 23. We first cut out 101x101 pixel 
postage stamps around each star In order to be sure the stars 
were centered, we computed the difference between the cen- 
troid of the star and the center of the array. If the difference 
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Fig. 7. — Measured curves-of-growth from aperture photometry for the objects in our sample in blue and red for F658N and F775W, respectively. The shaded 
regions show the 1 a uncertainties on the CoGs. The colored dashed lines denote the half-light radii, with the cross-hatched regions denoting the 1 a uncertainties 
on the radii. The colored dotted lines denote the resolution limit of the given image derived from the sizes of stars in the images, while the dotted curves show 
the CoGs of the image PSFs. The CoGs of Lya appear more extended than those of the rest-frame UV continuum, resulting in larger Lya half-light radii. 



was more than 0.2 pixels in either direction, we subsampled 
the image by a factor of 10, and shifted the star by one pixel 
for each tenth of a pixel it was offset from the center (the sub- 
sampling was done using the IDL function f rebin, which 
uses bilinear interpolation). The image was then binned back 
down to the native resolution. This process was run itera- 
tively on each star until they were all < 0.2 pixels from the 
array center Each star was normalized to its peak flux. The 
PSF was then calculated as the median of the five stars at each 
pixel position, and then normalized to a total flux of 1 . Mea- 
suring the half-light radii of the PSF in each band in the same 
manner as above, we measure an image resolution of r/, = 
0.09" in F658N, and r,, = 0.08" in F775W. Objects with r,, 
at or less than these values are considered unresolved at the 
limit of HST+ACS in their respective bands. 

For our sample of LAEs, we found half-light radii in the 
F658N image of 0.11" ± 0.02", 0.21" ± 0.02" and 0.20" 
± 0.02" for CHa-2, CH8-1 and CH8-2, respectively. In the 
F775W image, we found LAE half-light radii of 0.10" ± 
0.01", 0.16" ± 0.01" and 0.10" ± 0.01" for CHa-2, CH8- 
1 and CH8-2, respectively. Comparing these sizes to the PSFs 
discussed above, we find that we can definitively resolve 2/3 
LAEs in the F658N image, and 1/3 LAEs in the F775W im- 
age. At z = 4.4, the angular scale is ~ 6.671 kpc arcsec"' 
(for our assumed cosmology); thus these sizes correspond to 
0.7 - 1.4 kpc in Lya, and 0.7 - 1.1 kpc in the rest-frame UV 
continuum . 

4.2.2. Stacking Analysis 

As shown in Table 1, the signal-to-noise of the individual 
detections in the F658N image are not large. Thus, in order to 
obtain a more robust estimate of the average hal f-hght radii 
of LAEs, we have performed a stacking analysis. iHathi et aLl 
(|2008a) show in detail how such image stacking is justified 
for similar galaxies at similar redshifts, using the HUDF B, V 
and i' dropouts at z = 4, 5 and 6, respectively. 

Using the cutout stamps described in the above section, we 
first centered each LAE on the central pixel of each stamp us- 

^ The Year 7 Wilkinson Microwave Anisotropy Probe cosmology (Ho = 
70 km s"' Mpc"', Q,,, = 0.27 and Qa = 0.73 (Komatsu et al. 2010)) gives an 
angular scale of 6.899 kpc arcsec"', which would give physical sizes 4.3% 
larger than our assumed cosomology. 



ing the iterative technique described above for the PSF, requir- 
ing the SExtractor-derived center to be within 0.2 pixels of the 
center of the stamp. This step was performed separately for 
each object for each band, such that the F658N stamps were 
centered on the F658N emission, and the F775W stamps were 
centered on the F775W emission. Each centered LAE stamp 
was then normalized to its peak flux. A stacked image was 
then created in each band by taking the median of each pixel 
value from all three LAEs. Figure [8] shows the stacks of the 
three LAEs in both bands, with the contours denoting levels 
of constant brightness, as well as three-dimensional surface 
brightness profiles. 

We measured half-light radii of each of the two stacks (one 
for each band) in the same manner as the above section. The 
results from this analysis are shown in Figure |9] These stack- 
ing results confirm our observations of the individual objects, 
in that the Lya emission is more extended than the rest-frame 
UV continuum emission, with r/, = 0.16 ± 0.01" in F658N, 
and Th = 0.10 ± 0.01" in F775W. These angular sizes corre- 
spond to physical half-Hght radii of 1.07 ± 0.08 and 0.67 ± 
0.05 kpc for the F658N and F775W emission, respectively. 

5. DISCUSSION 
5.1. Rest-Frame UV Emission 

Inspecting Figure |7] one can see that LAEs are compact in 
their rest-frame UV continuum with half-light radii of r/, < 
1.1 kpc in all three objects. Thi s is consistent with pre vious 
studies of high-redshift galaxies. iFerguson et al.l (120041) stud- 
ied the rest-frame UV sizes of Lyman-break-selected galax- 
ies (LBGs) at z > 3, and photometric -redshift-selected galax- 
ies at 1 < z < 3 using data from HST. They found half-light 
radii from 0.25 - 0.4" at z > 2 (~ 2 - 3 kpc), rising to r/, ^ 
0.65" at z ~ 1 ('- 5 kpc) for galaxies with 0.7L* < Lj/y < 5L*. 
Similar size evolu tion has been fo und to extend out t o z ^ 6 
(IHathi et al.ll2008bi) and z = 7 - 8 dOesch et al.ll20Tol) . where 
Luv ^ L* LBGs have r;, ~ 1 kpc. 

Relativ ely few LAEs have had their morphologies studied. 
Recently, iBond et al . (2009) studied the rest-frame UV mor- 
phologi es of a sample of L AEs at z ~ 3.1 from the MUSYC 
survey dGawiser et al.ll20()6i) . They found that LAEs are typ- 
ically at least as compact as LBGs, with r/, < 2 kpc, and that 
the Lya emission is likely coincident with the UV emission 
(within < 1 kpc). .Gronwall et al.. (.2010.) studied the same 
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sample, examining the better detected LAEs (S/N > 30) in 
greater detail, finding that their rest-frame UV light is very 
concentrated, and that they have Sersic indices indicative of 
disk-like morphologies in most instances (0 < n < 2). 

5.2. Lya Emission 
5.2.1. Previous Results 

Prior to this study, only two high-redshift Lya-selected 
galaxies have been imaged in their Lya light at HST resolu- 
tion (i.e., using space -based narrowband da ta), pubhshed re- 
cently bvl Bond et al] (12010). In this studv. iBond et all (120 1 01 ) 
obtained HST/WFPC2 F502N imaging of z ~ 3.1 LAEs, ob- 
taining detections of two out of the eight LAEs they targeted. 
They concluded that these objects have Lya half-light radii 
< 1.5 kpc, similar to their rest-frame UV sizes, with the Lya 
emission coincident within 0.5 kpc of the rest-frame UV emis- 
sion, [khoads et al. (2009) also examined the relative sizes of 
LAEs in Lya and the UV continuum using ACS grism spec- 
troscopic data from the PEARS survey by examining the sizes 
of the objects in the spatial dimension. They did not find evi- 
dence of an extended Lya halo in a stack of the spectra from 
all 39 z ~ 5 galaxies in their sample. However, when stack- 
ing only the 10 galaxies with Lya observed in emission, they 
found that the spatial width of the spectrum at the position of 
Lya had FWHM = 0.26", while the same measurement on the 
adjacent UV continuum yielded FWHM = 0.19", suggesting 
possible extended Lya halos in these objects. 

Our positive detections of the three spectroscopically con- 
firmed LAEs more than doubles the total number of high- 
redshift LAEs with high-resolution imaging of their Lya light. 
Investigating the Lya light profiles of our LAEs, we find that 
the Lya emission appears relatively compact as well, with the 
half-light radius in every object at < 1 .4 kpc, and the mean 
size of '--^ 1.2 kpc consistent with th e Lya sizes of the two 
galaxies studied by IBond et al.l(l20I0l) . 

5.2.2. Individual Objects 

Comparing the CoGs of the Lya and rest-frame UV contin- 
uum light in individual LAEs in Figure |7] we find that CH8-1 
and CH8-2 have Lya half-light radii larger than the rest-frame 
UV, while CHa-2 is near the limit of our resolution in both the 
Lya and UV continuum light. The CoGs of CHS- 1 and CH8-2 
are very similar, with the rest-frame UV (F775W) profile ris- 
ing quickly, reaching the "total" flux at a radius smaller than 
the Lya (F658N) CoG, which is rising more slowly. Exam- 
ining the uncertainties on the CoGs, the difference between 
the Lya and rest-frame UV is at a > 1 sigma significance for 
much of the profile for CH8-2, and ^ I a for CH8-1. The un- 
certainties on If, are also quite small, with the half-light radii 
for Lya being 2 - 4 ct larger than that for the rest-frame UV 
for CHS- 1 and CHS-2. 

Investigating Figure |7] it is apparent that the low signifi- 
cance of the CHa-2 detection is hindering our measurement 
of its CoG, and thus its half-light radius measurement. Addi- 
tionally, for the remaining two objects, while their CoGs in- 
dicate larger half-light radii in Lya than in the rest-frame UV 
continuum, one will notice that their F65SN CoGs continue 
to increase out to the maximum radius. This effect is due to 
the CTE contribution to the background, which is a primarily 
positive signal caused by the overlapping CTE tails from the 
plentiful cosmic rays. It is thus possible that this CTE effect is 
artificially increasing the radii we measure in the F65SN data. 

It is thus prudent to examine these data to ensure that the 



result of larger sizes in the F658N data is a physical effect, 
and not an artifact of the data. We have performed a check 
on our results by measuring the sizes of galaxies that have 
F658N magnitudes similar to those in our sample, of 24.6 
< i^Fbsm < 25.2, yet have no Lya emission. For this test 
sample, we also required that the objects be detected at 5 u 
significance in both F658N and F775W, that 22 < nifj-i^w < 
29, and that the difference between the F658N and F606W 
magnitudes be < 0.1 mag. Out of our whole F658N-selected 
catalog, this yielded 67 objects. We further culled the sample 
by excluding objects near image edges, as well as highly ex- 
tended or clumpy objects, leaving a final sample of 28 objects. 
We measured the sizes of these objects in a similar manner as 
the LAEs in our main sample. We found the median of the ra- 
tio of Th.Fesm^^h.Fnsw to be 1.15, with a standard deviation of 
0.33. However, the uncertainty on the radii is much higher in 
the higher ratio objects; thus we computed a weighted mean, 
finding ih.F65m/i'h,Fn5w = 1.08 ± 0.02. This analysis shows 
that there is a slight systematic effect increasing the radii for 
objects in the F658N images over the F775W images. How- 
ever, with the exception of CHa-2 (which has the least sig- 
nificant detection, and thus is the most difficult to make con- 
clusions about), this ^ 10% effect is small when compared to 
the ratio of the radii for our three LAEs (1.1, 1.3 and 2.0 for 
CHa-2, CH8-1 and CH8-2, respectively) and the stack (1.6). 
We conclude that larger F658N sizes in our sample are likely 
real, but a larger sample of LAEs would increase the confi- 
dence in our result. 

5.2.3. Stacking Analysis 

As is shown in Figure [8] stacking the objects helps to re- 
duce the CTE-affected background. We see the same results 
in the stacking analysis in Figure|9]as hinted at in the individ- 
ual objects, with the CoGs of the rest-frame UV exhibiting a 
significantly steeper profile than that of Lya, highlighted here 
by the smaller uncertainties on the profile due to the greater 
signal-to-noise of the stacked images. Similar to the individ- 
ual results, the Lya half-light radius of the stack of LAEs is 
significantly greater than that of the rest-frame UV contin- 
uum, at ~ 5 <T significance. Also of note is that both CoGs 
reach a value of 1 by ^ 0.4", and oscillate around 1 (due to 
image noise) at higher radii. This implies that our derived 
half-light radii do not depend on our choice of a maximum ra- 
dius. We verified this, as changing the maximum radius from 
0.4 - 0.8" changed the resultant half-light radii by less than 1 
a. 

Figure|8]highlights this result, showing the two stacked im- 
ages with contours of constant brightness at 30, 50, 70 and 
90% of the peak flux. In the right-hand panels, we show 
three-dimensional surface brightness profiles of these images. 
As is evident to the eye, the contours on the Lya image 
are more loosely packed, and the Lya 3D image exhibits a 
broader slope than that of the rest-frame UV continuum im- 
age. Though the difference is slight, primarily due to the faint 
nature of these objects and the difficulty of space-based nar- 
rowband observations, these results are significant. 

5.2.4. Potential Future Improvements 

After the completion of our analysis, we were made aware 
of potential future imp rovements to the correctio n of the poor 
CTE in the ACS data dAnderson & Bedinll201(ih . While our 
analysis shows that the charge trailing is not significantly af- 
fecting our size measurements in the F658N data (§5.2.2), 
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Fig. 8. — Left: Stacked images of the three LAEs in the F658N band (top) 
and the F775W band (bottom), with both images normalized to their peak 
flux. The light in the F658N image is primarily due to Lya emission in these 
galaxies, while the F775W band shows the rest-frame UV continuum emis- 
sion longward of Lycj. The contours represent regions of constant flux, cor- 
responding to 0.3, 0.5, 0.7 and 0.9 of the peak flux. The F658N contours are 
also shown in the F775W image as dashed curves. Right: Three-dimensional 
surface profile of the F658N stack (top) and the F775W stack (bottom). The 
contours are the same as in the left-hand panels. The F775W stack has a 
steeper profile (and thus a smaller half-light radius), as shown by the more 
compact surface profile, and the denser contours. 

we were able to reprocess sections 13 and 22 of our F658N 
data using the updated CTE correction, in order to verify our 
results (see Figure |4|i. Briefly, the correction is based on a 
study of the trails behind warm pixels in dark exposures. The 
algorithm performs a mild deconvolution to restore the flux 
from the trails into the delta-function warm pixels. The cor- 
rection has been demonstrated to work well on backgrounds 
greater than 5 electrons, but at the time of development, suf- 
ficient data did not exist to calibrate the correction for back- 
grounds below this (the F658N images discussed here have 
backgrounds much less than this). Nevertheless, the current 
algorithm has been shown to correct the majority of CTE blur- 
ring, even at essentially zero background^. In an aperture of 
radius 0.3", the CTE-corrected data are ^ 0.2 mag deeper 
than the uncorrected data, which pushes these HST narrow- 
band data deeper than the existing ground-based data. 

For our three LAEs, we found F658N sizes in the CTE- 
corrected data of 0.18" ± 0.02", 0.21" ± 0.02" and 0.21" 
± 0.02" for CHa-2, CH8-1 and CH8-2, respectively. This 
implies that our size measurements for CH8-1 and CH8-2 are 
likely not adversely affected by the CTE problems, and also 
that CHa-2 may in fact be resolved, and larger in Lya than 
in the rest-frame UV. Stacking these objects, we find nearly 
identical results to our uncorrected stack, with r/, fdsm = 0. 17" 
± 0.01" and ih^fhsw = 0.10" ± 0.01". 

We currently plan to re-process all of our F658N data with 
the CTE correction, and do a new selection for LAEs based 
solely on the HST data to increase our sample of LAEs with 
resolved Lya emission. However, the CTE-correction algo- 
rithm needs to be verified at the low sky levels present in our 
data, and this requires new dark frames to be obtained. This 

^ One additional issue of CTE that enters in at low background is the im- 
pact of the read-noise, which did not go through the charge-transfer process. 
We examined the c orrections with an d without the readnoise mitigation em- 
ployed in Anders on & Bedinl <2010l) and found the resulting images to be 
essentially the same. 
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Fig. 9. — Measured fight profiles from aperture photometry for the stack 
of our sample of three spectroscopically confirmed LAEs in blue and red 
for F658N and F775W, respectively. The lines, shaded and cross-hatched 
regions are the same as in Figure|2] Both stacks exhibit the same result as the 
individual objects - Ly« emission appears more extended than the rest-frame 
UV continuum. 

work will be presented in a future paper. 

5.2.5. Interpretation of Results 

Our results indidicate that in our sample of LAEs, the Lya 
light is emitted from a larger region than the rest-frame UV 
continuum light. This result is intriguing, since both types of 
photons likely originate in the same location — the H II re- 
gions within the galaxy — thus one may expect both sets of 
photons to exhibit the same light profiles. However, Lya pho- 
tons are resonantly scattered by neutral hydrogen, while the 
rest-frame UV continuum is not. In an interstellar medium 
(ISM) that is homogeneous, there is no dust, this resonant 
scattering will result in a decoupling between the observed 
location of Lya emission and the rest-frame UV continuum 
emission, with much of the Lya emission eventually escap- 
ing from a random location far from its origin, appearing as 
an extended halo. However, in recent years, we have learned 
that many LAEs do in fact contain dust (e.g., [Pi rzkal et aj] 
l2007HFinkelstein et al.ll2008ll2009l:lPentericci et aLi200 9). In 
a dusty homogeneous ISM, where dust is evenly mixed with 
neutral hydrogen, resonant scattering will result in the major- 
ity of Lya photons being absorbed by dust; thus a pure ho- 
mogeneous ISM is unlikely, given that these galaxies exhibit 
Lya in emission. 

On the other hand, if the ISM is inhomogeneous, Lya 
can still escape in a scattered halo even if dust is present 
(Neufeld 1991; Hansen & Oh 2006; Finkels tein et"al] 120071 
2008, 2009). In an idealized case where the ISM is purely 
clumpy, nearly all Lya photons can escape as they scatter off 
of H I at the surface of the clumps and are thus screened from 
the dust. Even if an ISM is only partially clumpy, this ge- 
ometry still allows Lya to escape, though predominantly in 
a scattered halo. Additionally, if the majority of Lya emis- 
sion escapes only after scattering off of the receding edge of 
an outflow, the (now redshifted) Lya photons would also be 
decoupled from the rest-frame UV continuum, perhaps ap- 
pearing in a larger halo as well. 

Such halos have been predicted, but have yet to be observed 
at high-redshift. Given the modest signal-to-noise of our Lya 
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Fig. 10. — Filter profiles of the ground-based NB656 and NB665 filters 
are shown as the solid black lines, while the HST/ACS F658N filter profile 
is shown as the dashed black line. The position of the Lya emission fines 
for CHa-2, CH8-1 and CH8-2 are shown by the blue, red and green fines, 
respectively. The colored shaded regions denote the extent of the emission 
fines, as measured from the line FWHMs from the observed IMACS spectra. 

imaging detections, it is likely that we have not detected the 
full extent of these halos — rather we are just seeing the tip 
of the iceberg, in that the Lya light appears more extended 
that the rest-frame UV continuum. However, given the low 
signal-to-noise of the LAEs in our data, it is difficult to tell 
if we are truly seeing the edge of the Lya emission. One 
way to check whether our imaging has captured all of the Lya 
light is to compare the fluxes from the F658N imaging to that 
from ground-based photometry, which can be more sensitive 
to diffuse emission given the larger telescope apertures, and 
reduced sensitivity to read noise. 

We can perform this analysis for our sample of LAEs, 
which have their ground-based narr owband magnitudes tabu- 
lated in Table 1 of lFinkelstein et all(l2009i) . These magnitudes 
are 24.15 ± 0.11, 24.44 ± 0.15 and 24.39 ± 0.16 for CHa- 
2, CH8-1 and CH8-2, respectively. Comparing these mag- 
nitudes to those of the same objects from HST in the F658N 
data, we find that all objects have significantly greater ground- 
based narrowband fluxes than from HST, with flux ratios of 
iground/iHST of 2.0 ± 0.6, 2.0 ± 0.5 and 1.3 ± 0.3 for these 
three objects, respectively. This provides further evidence that 
these objects have significant Lya emission escaping in a dif- 
fuse halo, as the ground-based narrowband imaging detects 
2x more flux than HST for two objects, and 30% more for 
one object. However, this is not quite so straightforward, as 
the interpretation depends on exactly where the Lya flux re- 
sides in their respective narrowband filters. 

Figure [To] shows the ground based narrowband filters used 
to select these three objects, as well as that of F658N on ACS. 
The colored vertical lines denote the position of Lya for these 
three objects, while the lighter shaded regions denote the full- 
width at half-maximum of these lines, as measured from the 
IMACS spectra presented in §3.1. From this figure, we can 
see that all objects have their Lya emission lines encompassed 
within the FWHM of the F658N filter However, only CHa-2 
has Lya within the FWHM of the NB656 ground-based filter 
— both CH8-1 and CH8-2 were detected even though their 
redshift puts the emission line at < 30% of the peak F665N 
filter transmission. This means that the factor of two flux in- 
crease in CHa-2 from the ground over HST appears real, while 
the factor of 2 and 1.3 for the other two objects are actually 



lower limits, implying that the extended Lya halos in these 
objects may contain a dominant fraction of the total escaping 
Lya flux. 

We caution that as these objects are near the image depth 
limits in both sets of data, there could be zeropoint issues. As 
a test, we compared the narrowband fluxes from the ground- 
based NB665 ir nage (using 2.3"- aperture magnitudes from 
the catalog from lFinkelstein et alj (12008) ) to the F658N data 
to see if there is a zeropoint offset. We examined objects in 
common in both catalogs, computing the mean magnitude dif- 
ference in bins of 0.5 mag. From 21 < nifdsm < 25, the mean 
magnitude difference is always < 0.1 mag. However, there is 
of course scatter in the individual objects. At bright magni- 
tudes, this is small, with a ^ 0.3 mag, increasing to cr ^ 0.5 
mag at m^ sssAf < 25. The magnitude differences for 2/3 LAEs 
is significant, as their differences are greater than the 1 a un- 
certainty in the relative zeropoints. Thus, this flux difference 
is intriguing, and a similar analysis with a larger LAE sample 
will provide greater confidence in this effect. 

6. CONCLUSIONS 

We have performed high-resolution HST F658N narrow- 
band imaging over a portion of the GOODS CDF-S in order to 
directly image resolved Lya emission at z = 4.4. We have de- 
tected Lya emission from three spectroscopically confirmed 
LAEs in these data. 

Studying the relative positions of these objects in their Lya 
light (from the F658N data) and their rest-frame UV contin- 
uum light (from existing F775W data), we find that our data 
do not support a positional offset between the two types of 
emission. We then measured the light profiles and half-light 
radii from our three LAEs in both filters. We find that in all 
three objects the Lya light profile rises more slowly, and has 
a larger half-light radius than the rest-frame UV continuum 
emission. We confirmed this result by stacking the galaxy im- 
ages in both bands, finding that the Lya emission has r/, = 
1 . 1 kpc, while the rest-frame UV continuum is more compact 
with r/, = 0.7 kpc. This implies that the Lya light is more 
spread out, presumably due to effects of resonant scattering, 
possibly in a clumpy ISM. 

To investigate this further, we compared the fluxes of our 
LAEs in the F658N narrowband to ground-based narrowband 
measurements, which are more sensitive to diffuse emission 
due to, among other things, larger telescope apertures. We 
find that in all three cases where we have measurements from 
both ground and space, the ground-based narrowband fluxes 
are significantly greater than the space-based fluxes, by fac- 
tors of 1.3 - 2.0. This shows that the larger physical sizes 
detected in the F658N data are only the tip of the Lya ice- 
berg, and that the majority of the Lya emission may lie in a 
larger, diffuse halo. It is thus clearly important to include the 
ISM geometry and kinematics in any study of Lya emission 
at high redshift. 

While Lya emission is one of the most powerful tools we 
have to discover and study galaxies at high redshift, the com- 
plicated radiative transfer undergone by Lya photons in their 
host galaxies muddle the physics that can be inferred. In or- 
der to maximize our understanding of LAEs and Lya emis- 
sion in general, we need to obtain a greater understanding of 
how Lya makes its way from the H II regions where it origi- 
nates to its point of escape from the galaxy. Studying the Lya 
spatial profiles provides one estimate of the complex radiative 
transfer by comparing the Lya morphologies to those of the 
rest-frame UV, but more work is needed to obtain strong de- 
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tections of these diffuse halos, which likely requires the next 
generation of ground and space-based observatories. 
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